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Abstract 



The thermoelectric power S and the thermal conductivity k of stoichiomet- 
ric CaBg and vacancy-doped Cai-^Bg have been measured between 5 and 
300 K. The thermopower of both materials is surprisingly large at room tem- 
perature. Across the whole covered temperature range, S is negative and 
the temperature dependence is most likely dictated by band-structure effects. 
The empirical interpretation of our data involves a calculation of the ther- 
moelectric power by using the Boltzmann equation in the relaxation-time 
approximation and the assumption of a defect band in the proximity of the 
lower edge of the conduction band. Good agreement with our data is found 
by considering acoustic-phonon and ionized-impurity scattering for the elec- 
trons in the conduction band which is well separated from the valence band. 
The thermal conductivity k{T), predominantly due to phonons, may also be 
described quite well across the whole covered temperature regime on the basis 
of a Debye-type relaxation-time approximation and assuming the concurring 
influence of boundary-, dislocation-, Rayleigh-, and resonant-scattering on the 
mean free path of the phonons. 
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I. INTRODUCTION 



Recent experiments on hexaborides with divalent metal cations of the alkaline-earth 
series, namely Ca, Sr and Ba, have revealed some unusual physical properties of these mate- 
rials. In particular, an itinerant type of weak ferromagnetic order, stable up to temperatures 
between 600 and 900 K, has been observed in CaBg and related alloys in a very narrow range 
of electron doping.&i The electronic properties of alkaline-earth hexaborides place these ma- 
terials close to a metal-insulator transition. ii This has been confirmed theoretically by recent 
calculations of the electronic structure of divalent hexaborides which have shown that, with 
the exception of a small region around the X point of the cubic Brillouin zone, the valence- 
and the conduction bands are separated by a large gap of several eV.i The slight band 
overlap is reduced if the inter-octahedron B-B bond distances increase, finally resulting in 
the opening of a gap over the entire momentum space. It has also been predicted that the 
transport properties of SrBg may be strongly dependent on doping, with changes from a p- to 
an n-type metal induced by slight shifts of the Fermi energy around the zero doping value.ll 
These observations have led to the speculation that the ground state of undoped divalent 
hexaborides may be characterized by a Bose condensate of bound electron-hole pairs, or 
excitons.iHi Weak ferromagnetism may then develop because of a spontaneous time-reversal 
symmetry breaking via doping. More recent calculations of the single-particle excitation 
spectrum of CaBg, however, have led to the claim that this material is not a semimetal 
but a semiconductor with a minimum band gap of 0.8 eV.i If this possibility is consid- 
ered, the experimental observations, which suggest that binary hexaborides are close to a 
metal-insulator transition, would indicate the presence of a defect band of itinerant charge 
carriers. In an attempt to resolve the uncertainties with respect to the correct description 
of this material, we have studied the influence of the chemical composition on the physical 
properties of CaBg. The electrical resistivity, magnetoresistance, low-temperature specific 
heat and the optical conductivity of stoichiometric and doped CaBg samples have been 
studied and the results have been presented in a previous publication.! Below, we present 
the results of measurements of the thermal conductivity and the thermopower. While data 
on the electrical conductivity and the electronic contribution to the thermal conductivity 
provide information on the density of electronic states N{E) close to the Fermi energy Epi 
the thermopower may serve to establish the energy derivative of N[E), i.e., an additional 
important detail of the electronic excitation spectrum. One of the samples was prepared 
in such a way as to obtain material with a close to stoichiometric composition, denoted 
as CaBg. The second sample, which we denote as Cai.^Bg, contained a small number of 
vacancies on the calcium sites, inadvertently introduced during the flux-growth procedure 
and leading to a certain degree of self doping which is difficult to control. This paper is or- 
ganized as follows. After a brief description of the sample preparation and the experimental 
methods used in this investigation in section we present, in section the results of our 
measurements and their analysis. The conclusions are presented in section 



II. SAMPLES AND EXPERIMENTAL METHODS 

Binary hexaborides can be synthesized in a narrow range of composition, with the 
tendency to be boron-rich with concentrations of metal vacancies up to several percent. 
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Stoichiometric CaBg crystals can be obtained close to the border of the metal-rich phase 
boundary.!'! For our experiments the single crystals were grown by a slow-cooling procedure 
in aluminum fluxlll starting with a nominal ratio CaBs for the stoichiometric CaBg sam- 
ple, and with CaBi2 for the metal- deficient Cai.^Bg sample. The crystals were removed 
from the flux by leaching in a concentrated sodium hydroxide solution. Subsequent etching 
with HNO3 was intended to remove possible surface contaminations. The samples on which 
our transport measurements were made were of prism-type shape with approximate overall 
dimensions of 4.2 x 0.5 x 0.45 mm^ for CaBg and 4.7 x 0.45 x 0.4 mm^ for Cai.^Bg. 

The thermoelectric power S and the thermal conductivity n of both samples were mea- 
sured simultaneously by means of a standard steady-state heat-fiow technique. A commercial 
^He gas-fiow cryostat was used for cooling the sample holder. At one end of a prism-shaped 
sample, the thermal contact to a copper heat sink was achieved by using high- conductance 
silver epoxy. The sample heater, consisting of a 100 Q ruthenium-oxide chip resistor, was 
attached to the other end of the prism by the same method. Joule heating caused by heater 
currents of the order of a few mA provided the necessary heat-flow and hence a thermal gra- 
dient along the crystals. In order to measure both the temperature difference AT between 
two contacts mounted perpendicularly to the heat flow as well as the thermoelectric voltage 
of our samples with respect to chromel, we used two pairs of calibrated 0.025 mm Au-Fe 
(0.07 at. %) versus chromel thermocouples.0 In order to minimize spurious thermal volt- 
ages, the thermocouple leads were connected to uninterrupted copper wires reaching three 
home-built low-noise voltage amplifiers, mounted directly on top of the cryostat insert. 



III. EXPERIMENTAL RESULTS AND ANALYSIS 



A. Thermoelectric power 

The thermoelectric power S{T) of both samples is shown in Fig. ^on linear scales. The 
negative sign of S clearly demonstrates the dominant n-type character of the investigated 
crystals. The rather large values of S, if compared to the typical values observed for the 
thermoelectric power of common metals, are an indication for the low itinerant charge- 
carrier concentration in these materials, compatible with the very low values of the electrical 
conductivity reported previously.i The overall temperature dependence of S is obviously non- 
linear, with features that are reminiscent of band-structure effects affecting the thermopower 
of elements and crystalline alloys.^ Nevertheless, below T* = 20 and 40 K for CaBg and 
Cai-^Bg, respectively, the curves show an approximately linear temperature dependence. A 
negative thermoelectric power varying linearly with T is usually identified as the diffusion 
thermopower 5* of metals, where the Fermi-energy Ep is located within the conduction band. 
In the free-electron approximation, 

^~ 3\e\Ep/s' 

with E-p as the Fermi-energy measured from the bottom of the conduction band, and s a 
factor, typically of the order of unity, describing the energy dependence of the scattering 
time r oc E^~'^/'^ for the scattering mechanism dominating the low temperature behaviour .llHl 
From the average slope of S{T) below T*, we may thus calculate a rescaled Fermi energy 
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E-p/s. These values are E-p/s = 4.5 and 8.7 meV for CaBg and Cai-^Bg, respectively. 
Considering the entire temperature dependence of S, in particular the abrupt changes of 
slopes of S{T) around T*, suggests that ultimately, S(T) cannot be explained by simply 
taking into account electronic states in the conduction band alone. In the following, we try 
to identify the possible reasons for the measured departure from linearity in the temperature 
dependence of S. 

With the common assumption that the current densities j and w for electrical and 
thermal transport, respectively, respond linearly to their driving forces, we can write 

j = e%E-^LiVrT, 

1 (2) 
w = eLiE - -LsVrT, 

where E is the external electric field and VrT is the temperature gradient. In Eq. H, the 
transport coefficients Li (i=0,l,2) explicitly fulfill the Onsager relations.^ For j = 0, we 
obtain w = — KdVrT and E = SV^T, whereas for VrT = 0, we get j = crE. These 
are the equations defining the electronic contribution to the thermal conductivity Ke\, the 
thermoelectric power S, and the electrical conductivity a, such that 



1 / L^' 

Kei = ^ f ^2 - ) , (3a) 



\e\T Lq 

a = c^Lq. (3c) 



The integrals Lj (i=0,l,2), are defined as 

U = - d{E){E-Epf-^dE, (4) 

with a{E) representing the conductivity spectrum, including all system-dependent features. 
By partial integration of Eq. and by using the Sommerfeld expansion scheme, we get 



3|e| V dE 



1 



(5) 



which approximates Eq. |3b| very well at temperatures ksT -C Ep. Since we intend to 
interpret our data over the entire covered temperature range, we will not use the low- 
temperature approximation represented by Eq. H, in our analysis. For our purposes the 
general form of S shown in Eq. ^ is more suitable. 

Since o-{E) > 0, the coefficient Lq is strictly positive and the sign of S is determined by 
the sign of the integral Li. By inspecting the integrand of Li, represented in Eq. H, and 



also by considering Eq. SB, it may be concluded that the states with energies higher than 
the chemical potential, i.e., E > Ep, provide a negative contribution to the thermoelectric 
power, whereas those states which are located below Ep contribute with the opposite sign. 
At low temperatures, where only the states very close to Ep contribute to the electronic 
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transport, the measured 5* values are negative. This leads us to conclude that d'{E) must 
increase monotonously across Ep, so that the states contributing negatively to S(T) in 
Li acquire a larger weight than those located below the Fermi-energy, which contribute 
positively to S{T). In order to reproduce the observed reduction of the slope of S(T) at T* 
in the calculation, however, we are forced to introduce a "feature" in cr(i?) centered at an 
energy E* such that \E* — i^pl ~ kBT*. Since the absolute value of dS/dT decreases above 
T*, the feature in a{E) must either reduce the negative contribution to S{T) given by the 
states above E-p-, or else enhance the positive contribution to S{T) of the states below E^- 
In what follows, we assume that in our case the second possibility is realized and we add an 
additional part to the cf{E) spectrum, decreasing linearly with increasing energy below E^- 
The complete chosen (j{E) spectrum, corresponding to the Eqs. |^ and |10| below, is shown 
in Fig. |. 

Based on this preliminary analysis and the model outlined above, we have calculated 
the electrical transport properties of CaBg and Cai-^Bg. The Fermi energy is fixed within 
a parabolic conduction band. Close to the bottom of this band we position an additional 
band, which may be interpreted as a defect band. With these assumptions, the conductivity 
spectrum may be written as 

d{E) = d^{E) + dP{E). (6) 

For an isotropic parabolic conduction band, the Boltzmann equation in the relaxation 
time approximation leads to 

where the bottom of the conduction band has been fixed to = meV. By setting 
cf{E) = a^{E) in Eq. ^, it is easvto recover Eq. |1|. For the effective mass of the conduction 
electrons we used m„ = 0.28mopi with tuq as the free electron mass. We only considered 
two scattering mechanisms of the charge carriers in the conduction band, i.e., scattering by 
acoustic lattice vibrations, approximated by a rateHll 



r, 



Iv 



aE^/^T, (8) 



and scattering by ionized impurities, with a rateli3 

r,T' = PE-'/\ (9) 

where a and (3 are two constants to be determined by the fitting procedure. If the defect 
states which are forming a band below the bottom of the conduction band coexist with 
localized states, the latter may act, when ionized, as scattering centers for the electrons. 
This justifies considering a term r^-^ in the scattering rate of the electrons. As we will 
show below, the existence of localized donor states may also be compatible with our thermal 
conductivity data. Another possibility is to attribute the rate t-^^ to the vacancies at the 
Ca^~'"-ion sites which would act as negatively charged scattering centers. With the common 
assumption that the two scattering mechanisms do not interfere with each other, we may 
use Matthiessen's rule in the form r = (t^^ + r^'^)^^ for the electrons' average scattering 
time. 
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For the defect band, we postulated 

\0 {E>E„,), ^ ' 

with a mobihty edge E^^ and 7 a constant. Since the quantities Lq and Li which dictate 
the temperature dependence of the thermopower, appear in the ratio Li/Lq, S depends 
only on the relative magnitude of the constants a, f3, and 7. This allows some reduction 
of the number of fitting parameters. In Fig. 0, the solid curves represent the results of the 
calculation of S{T). By considering the simplicity of the model and the small number of 
free parameters (a/7, /?/7, Ep, and -Eme), the agreement with the experimental data is quite 
remarkable. The Fermi energies which emerge from the fitting procedure amount to 14 and 
28 meV for CaBg and Cai-^Bg, respectively, corresponding to K charge- carrier densities 
ric of 1.1 X 10^^ m~^ and 3.2 x lO^'' m~^, or to 0.8 x 10"'^ and 2.3 x 10""^ charge carriers 
per unit cell, respectively. Preliminary results of low-temperature Hall effect measurements 
performed on a CaBg sample with the same resistivity as our vacancy-doped sample lead, 
when interpreted with a single-band model, to a charge carrier density of 4.2 x 10^^ m~^ 
(Ref. |1^), hence compatible with our results. The value of the mobility edge -Eme, introduced 
to reproduce the kink in S{T) at the observed temperatures is approximately +8 meV for 
both samples. We note, however, that this value strongly depends on the choice of the 
energy exponent of the low-temperature scattering term. In this sense we have to admit 
that this result depends on the employed model. 

Our calculation of the thermoelectric power relies essentially on the model of the conduc- 
tivity spectrum (modulo a factor), which we may now test by calculating other electronic 
transport properties such as the electrical resistivity p(T) and the thermal conductivity 
/€ei(T) and compare the results with the experimental data. In Fig. ^, we have plotted the 
rescaled electrical resistivity of CaBg calculated with Eq. ^ together with the experimental 
data taken on two samples of the same batch. The agreement between the calculated curve 
and the experimental data is acceptable, especially if one considers primarily the tempera- 
ture dependence. 

By using Eqs. ^ and 0, combined with the measured resistivity curves shown in 
Fig. ^, we conclude that the mean free path Zp = hk-pT^E-p) /m^ of the conduction electrons 
in CaBg is less than 100 A. The same length is obtained, obviously enough, by employing 
the Drude formula a = ne^r/mn and by using n = Uc, the value cited above for the charge 
carrier concentration of CaBg. Since the electrons' wave vector at the Fermi-energy is of the 
order of 10^^ A~^, it seems fair to assume that the loffe's criterion for localization [kplp < 1) 
is close to be fulfilled in CaBg. 



B. Thermal transport 

In Fig. H, we show the measured thermal conductivities ^mcas of CaBg and Cai-^Bg 
on logarithmic scales between 6 and 300 K. In the further analysis we assume that only 
conduction-electron and lattice contributions to Kmeas need to be considered. 

A measure for the electronic contribution to the thermal conductivity is given by the 
quantity L* = K^i/aT, which can be calculated from the Eqs. ^ and Just like the 
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thermoelectric power 5", L* depends only on the relative magnitude of the constants a, /3, 
and 7. In a degenerate one-band model, the value 

nJaT = Lo = ^(^^y (11) 

is denoted as the Lorenz number. The ratio L* can be much higher in semiconductors 
than in metals, due to the so-called ambipolar diffusion, i.e., the diffusion of electron-hole 
pairs. In Fig. ^ we plotted the temperature dependence of L*, calculated by using the 
fitting parameters provided by our analysis of S{T), together with the value Lq = 2.45 x 
10^^(V/K)^. As expected, Kei/crT is equal to the Lorenz number at very low temperatures. 
This is so because only the electrons of the conduction band contribute to the transport. 
The influence of the ambipolar diffusion to the electronic thermal conductivities of CaBg 
and Cai-^Bg is manifest at intermediate temperatures. Nevertheless, this effect is still quite 
small if compared to enhancements of L* by two orders of magnitude observed in intrinsic 
semiconductors.lii 

The electronic contributions to the thermal conductivities of CaBg and Cai-^Bg, shown 
in Fig. ^ as solid lines, have beeen calculated by using the Wiedemann-Franz relation 
Kei = L*aT, with L* calculated from the Eqs. ^ and ^ as shown in Fig. ^. For the 
electrical conductivity a, we used the values which have typically been measured on close- 
to-stoichiometric and vacancy-doped samples, respectively.i'Eli^ For both hexaborides, the 
estimated Kei is at least two orders of magnitude smaller than /tmeas at all temperatures cov- 
ered in this study. Since the experimental uncertainty in our measurements of the thermal 
conductance is of the order of 0.5%, the procedure for estimating k^i is not critical and may 
be considered as sufficiently accurate for the validity of the following considerations. 

The lattice contributions /tph = ^meas — ^ei to the thermal conductivities of both samples 
are shown in Fig. ^ revealing monotonic increases with T up to 18 K and 23 K for CaBg 
and Cai_<5B6, respectively. At these temperatures, /tph passes over maxima of 124 and 74 
Wm~^K~^, respectively, and decreases with varying slopes up to room temperature. We 
note that at low temperatures Kph of Cai.^Bg attains lower values than /tph for CaBg. This 
seems reasonable because of an expected enhanced scattering of the phonons by defects in 
the Ca-deficient material. The crossing of the curves at a temperature close to 30 K is a 
rather surprising experimental observation which indicates a more efficient high-frequency 
scattering mechanism of the phonons in the structurally more perfect sample. 

The following quantitative analysis of our data is based on a relaxation-time approxima- 
tion in the context of the Debye model. The Debye approximation is valid at temperatures 
which are low enough with respect to the Debye temperature Gd, where we may assume 
that the acoustic branches uj{k) of the vibrational modes are well approximated by a linear 
k dependence. Since the temperature range covered in our experiment is considerably lower 
than the Debye temperature Bd = 783 K,i this approximation is justified to some extent. 
With this simplification, the thermal conductivity may be calculated from 



= 3nkB (^Q- J fph ^g^. _ ^y ^i^, T)dx, (12) 

where t{x,T) is the mean time between collisions of a phonon, n the number density 
of atoms in the crystal, and x = huj/k^T. By limiting the total density of allowed degrees 
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of freedom to 3n, it is possible to calculate the cut-off frequency ujd, which is related to 
the slope of uj{k). This cut-off frequency defines the Debye temperature Gd of the material 
which is, in this sense, a measure of the phonon velocity and is usually extracted from low- 
temperature specific-heat or elastic-constant measurements. For the Debye temperature 
6d, we used the above cited value which was extracted from specific- heat measurements 
on stoichiometric CaBg.i Since upon introducing a very small number of calcium vacancies, 
no significant changes in the lattice contribution to the specific heat are expected, we used 
the same Debye temperature for CaBg and Cai.^Bg. This argument is supported by the 
observation that there are no significant differences between the low-temperature specific 
heats of CaBg and Cai_^La2,B6 with x=0. 005.1 The number density of atoms in CaBg can 
be calculated by using the lattice constant a = 4.146 Areported in Ref. ^. From the relation 

n = l/67r2(A;BeD/Hh)', (13) 

we may also calculate the average velocity of the vibrational modes for which we obtain f ph = 
5722 m/sec. This value is compatible with the average sound velocity f ph = 4200 m/ sec, valid 
for EuBgii when scaled with the mass-density factor v^PeuB^7a>B6 • 

For evaluating the integral in Eq. [1^, we approximate the phonon relaxation rate by 



T-\x, T) = {t^I + r,Ti + r^l^^ + r'Dix, T). (14) 

The different terms on the right hand side of Eq. |l^ represent the scattering of phonons 
at grain boundaries, by the strain fields surrounding dislocations, by point defects and via 
resonant scattering, respectively. The resulting fits are shown in Fig. ^ as solid lines. In 
what follows, we discuss each term appearing in Eq. separately. The values of the free 
parameters involved in our fitting procedure (see below) are given in Table | and the resulting 
total scattering rates of the phonons are shown as a function of frequency in Fig. ^for T = 10 
and 100 K. 

The temperature and frequency independent scattering length /cas in the first term in 
Eq. |14|, TqI^ = fph/^cas, is the so-called Casimir lenght, where the phonon mean free path 
is limited by the sample dimensions. ill The values of 0.303 mm and 0.083 mm for CaBg 
and Cai-^Bg, respectively, are of the order of magnitude of the smallest sample dimensions, 
resulting from the variation of the sample cross section perpendicular to the direction of the 
thermal current. 

The second term in Eq. |l^, Tj^^j = A-cu, is usually attributed to phonon-electron scattering 
but, as pointed out in Ref. can also be caused by phonon scattering due to strain fields 
surrounding dislocations. As may be seen in Table |, this term is slightly larger for CaBg 
than for Cai-^Bg. In view of the higher density of charge carriers in Cai-^Bg derived in 
section p.Il7\|, this rules out that mainly phonon-electron scattering is responsible for this 



term. Phonon-electron scattering may also be excluded by inspecting the relevant scattering 
rate which can be written aS 



r. 



(15) 



with Cd ~ ~|-^F as the deformation potential. The fit values for A given in Table | result in 



Fermi-energies of the order of 1-2 eV, i.e., values which exceed Ep evaluated in section [III A 
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by two orders of magnitude. For phonon scattering at dislocations we haveE3 

r^, = 0.06nd(67G)'^, (16) 

with rid as the dislocation density, 7g the Griineisen parameter, and b the Burger vector. By 
using the lattice constant for 6 ~ 4 Aand 7g ~ 1, we obtain nd ~ 2 x 10^^ m~^, a plausible 
result. 

The enhanced presence of scattering centers at the vacancy sites in Cai-^Bg is unequiv- 
ocally reflected in the third term of Eq. |14|. Rayleigh scattering by point-like mass defects 
leads, in the continuum approximation, to a scattering rate term Tj^^yi ~ ^ ' ^'^^ 

with0 

In Eq. Hp is the number density of point defects in Cai-^Bg, a = 4.146 Ais its lattice 
constant, M is the mass of a unit cell, and 6M is the supplementary or missing mass due to a 
point defect in a sphere with the radius of the point defect itself. Since the extra mass term 
appears squared in Eq. |1^ it is not possible to decide whether 6M is positive or negative. 
By assuming that each point defect in Cai^s^e provides one free electron to the ensemble of 
conduction electrons and that for stoichiometric CaBg = 0, we have rip = 3.2 x 10^^ m~^, 
i.e., the same value as the electron density calculated in section [111 A|. Again, by using 



for B the value which we obtained from the fitting procedure, we note that SM/M is of 
order unity. In view of the rather crude approximations, this result is encouraging, as it 
seems to indicate that the concentration of itinerant electrons in CaBg is equivalent to the 
concentration of impurities (or vacancies) with approximately the same mass as the mass 
of the unit cells. The scattering rate term Xj^^yi is only of significance for fitting the fCph 
data of the vacancy-doped Cai.^Bg sample. Our attempts to evaluate the free parameter 
B for CaBg resulted in a vanishing contribution of r^lyi to the total scattering rate of the 
stoichiometric hexaboride. Considering the estimated error range of a few times 10~^^ sec^ 
for B, this is at least consistent with a reduced amount of defects in CaBg. 

The last term appearing in Eq. |T3| needs a detailed examination. Resonant scattering 
of phonons has already been invoked in the analysis of the thermal conductivity data of 
EuBg in a previous invest igat ion. il The appearance of in the scattering rate of EuBg 
was intepreted as being due to a resonant scattering of low-momentum vibrational states in 
the dispersionless part of the uj{k) branches with acoustic phonons of constant energy that 
seem to occupy an extended region of the Brillouin zone at higher momenta. Indeed, a rapid 
flattening of the LA branches at an energy of approximately 12 meV has been observed in 
inelastic neutron scattering experiments performed on the trivalent hexaborides XBg (X = 
Ce, La, Sm).&^i Because of similar atypical temperature dependencies of Kph above 50 K, 
we follow the arguments given in Ref. and introduce a resonant scattering term of the 



form 



^r-el=C-— -^F{T), (18) 

with 



^^^^ expihuj^JkBT) - 1 ^^^^ 
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as the bosonic population factor. In our intepretation of the data, the crossing of the thermal 
conductivity curves at T = 30 K has to be traced back to a much stronger resonant scattering 
rate in CaBg. As may be seen in Table |, the fit value C is more than one order of magnitude 
larger for CaBg than for Cai_^Bg. 

Although, by taking into account the resonant-scattering term, we achieve a fairly ac- 
curate description of our thermal conductivity data, it seems rather unlikely that the large 
difference in the temperature dependence of Kph may be accounted for by considering phonon- 
phonon scattering processes alone, especially if we consider the only small difference in 
chemical composition between the two materials. 

In order to find a more plausible explanation for the experimentally observed crossing of 
/tph(r), we note that a strong reduction of the lattice thermal conductivity of semiconductors 
induced by defect donor states has been reported in a large number of publications.!!^ The 
intensity of the phonon scattering processes invoking such states depends strongly on the 
electronic excitation spectra induced by these defects. In particular, donor states in cubic 
semiconductors, leading to the existence of localized singlet-triplet centers just below the 
conduction band, are also expected to provide the source for a strong resonant phonon 
scattering.i§ Physical realizations of such systems have been found in the valley-orbit-split 
ground state of the donors As, Sb, P, and Li in Ge. The singlet-triplet valley-orbit splitting 
of the fourfold degenerate ground state of As in Ge, for example, is 4.23 meV, a value which 
is of the same order of magnitude as the resonant energies of 10.85 and 9.38 meV, quoted 
in Table | for CaBg and Cai-^Bg, respectively. If the density of defects in the crystal is 
such that the average distance between the donor states is of the order of the extension 
of the defect wave function, characterized by a Bohr radius of typically 40 A, we expect a 
broadening of the degenerate energy-eigenvalues, finally resulting in the formation of a band 
of defect states. The typical extension of the defect wave function should be compared with 
the expected mean distance between the electrons in CaBg, estimated by calculating the 
third root of the inverse of the electron density derived in section [111 A| . This leads to values 
of 70 A and 100 A for Cai-^Bg and CaBg, respectively. If the electron densities at K 
derived in section III A are due to the existence of donor states in CaBg, we are thus close 
to concentrations for which a band is expected to form. A schematic representation of this 
is shown in Fig. Because of the statistical distribution of distances between donor states, 
it seems possible to realize a situation where both, localized donor states and defect-band 
itinerant states, coexist. Since stoichiometric CaBg is closer to perfection than the vacancy- 
doped sample, the observation of a stronger high-frequency scattering mechanism in CaBg, 
as mentioned above, might appear as rather surprising. However, as we noted above, an 
increasing distance between the defects due to a lower concentration of defects would narrow 
the tight-binding bands in CaBg, eventually resulting in a larger fraction of localized donor 
states acting as resonant scattering centers for the phonons. Hence, if this scheme is valid, 
the higher resonant scattering rate in CaBg makes sense. 

All our attempts to cast the scattering of phonons at localized donor states into a reason- 
able equation for its relaxation rate and to obtain a fair description of the high temperature 
features of Kph for CaBg failed. This failure may either be due to a breakdown of the Debye 
model at higher frequencies, or may indicate that localized donor states in CaBg do not 
scatter phonons. With respect to the first possibility it may well be that using the Debye 
model at the frequencies around the resonance is not justified. A strong resonant scattering 
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mechanism may, in some cases (see, for example, Ref. ^0]), modify the unperturbed disper- 
sion relations by opening a gap in the phonon spectrum at the resonant frequency, thus 
turning the Debye scheme into a very poor approximation at the frequencies around the 
resonance. This may also be the reason for the poor quality of the fit to the CaBg data, for 
which the resonant-scattering is quite significant. Since the dispersion relations of the lattice 
vibrational modes in CaBg have not been measured yet, we cannot offer a more detailed 
discussion. In case of deviations of the acoustic branches of CaBg from a non-dispersive k 
dependence at the frequencies which are relevant for our investigation, the functional form 
of r~g would certainly change. 

Concluding this section, we recall that our thermal conductivity data show a strong and 
somewhat unexpected dependence on calcium vacancy-doping at elevated temperatures, 
which we interpret as being due to a resonant scattering of the vibrational modes. For 
a more realistic computation of the thermal transport properties at the high end of the 
temperature regime covered in our investigation, more information concerning the phonon 
dispersion relations of CaBg is necessary, however. 



IV. SUMMARY AND CONCLUSION 



We have measured the thermoelectric power and the thermal conductivity of close to 
stoichiometric and vacancy-doped CaBg between 5 and 300 K. The high negative values of 
the thermoelectric power S{T) of both materials indicate a low concentration of itinerant 
n-type charge-carriers. We have achieved a reasonable interpretation of S{T) by using a 
relaxation-time approximation of Boltzmann's equation, by considering a sizeable gap be- 
tween valence and conduction band states, and by postulating the existence of an additional 
term in the conductivity spectrum which may be related to the existence of a defect band in 
proximity to the lower conduction band edge. A Debye-type relaxation-time approximation 
is successful in describing the thermal conductivity data across the entire temperature range 
covered in the experiments. The usual terms contributing to the total scattering rate of the 
vibrational states agree well with our expectations. A rather unusual temperature depen- 
dence of Kph is observed above approximately 30 K. A resonant type of scattering seems to 
influence the mean free path of the phonons at these temperatures. Quite surprisingly, its 
limiting effect on the time period between two vibrational-mode collisions is much stronger 
in the material closer to stoichiometry. Although postulating a resonant-type of phonon- 
phonon scattering leads to a satisfactory description of our experimental observation, an 
alternative possible interpretation which considers the existence of localized donor states 
is suggested. Since the lattice dispersion relations of CaBg have not yet been investigated 
experimentally, we cannot offer any rigorous statements about the physical nature of this 
resonant type of scattering of the vibrational modes. 
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TABLE I. Parameters of the fitting of the Kph(T)-data to Eqs. 



12 



and 



14. 



Fit parameter Units CaBg Cai-^Bg 

kjZ Jmi 303 ±4 83 ± 1 

A 10-6 20.5 ±0.1 18.0 ±0.5 

B 10-^"^ sec^ - 7.7 ±0.1 

C 10^0 sec-i 98.5 ±0.1 7.9 ±0.1 

niores meV 10.847 ±0.003 9.38 ± 0.02 
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FIG. 1. Temperature dependence of the thermoelectric power 5* of CaBe and Cai-^Bg on 
Unear scales. The solid lines represent fits to our data, as explained in the text. 
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FIG. 2. Auxiliary quantities, relevant in our analysis of S{T). The solid line is the postulated 
conductivity spectrum ^{E). The various broken lines represent the value of the integrand appear- 
ing in Li as a function of energy for T = 50, 100, 200 K, as explained in the text. In the inset, the 
integrand of Li is shown as a function of T and E. For a given temperature, the thermoelectric 
power S is proportional to the integral taken along the curves shown on the three-dimensional 
surface. 
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FIG. 3. The solid line represents the expected temperature dependence of the electrical resistiv- 
ity of CaBg which is based on the postulated conductivity spectrum a{E), as explained in the text. 
Also shown in the figure as dotted and broken line are the measured temperature dependencies of 
the resistivities of two samples originating from the same batch as the CaBe sample measured in 
this investigation. 
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FIG. 5. The ratio Ke\/aT as a function of temperature, calculated by using the Eqs. 3a and 
|3^ , with the fit parameters provided by our analysis of S{T) of CaBg (solid line) and Cai-^Bg 
(dashed line). The horizontal dashed-dotted line is the Lorenz number Lq = 2.45 x 10~^ (V/K)^. 
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FIG. 6. The lattice contribution to the thermal conductivity Kph of CaBg and Cai-^Bg as 
a function of temperature on logarithmic scales (on linear scales in the inset). The solid lines 
correspond to fits to the data as explained in the text. 
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FIG. 7. The total phonon relaxation rate r i resulting from our fit (see Eq. ^) as a function 
of frequency for T = 10 and 100 K. 
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FIG. 8. Schematic representation of the bandwidth dependence on defect spacing. A higher 
defect concentration in the vacancy-doped sample Cai-^Bg decreases the spacing between donor 
states, resulting in a higher concentration of overlapping states. The non-overlapping, localized 
states are postulated to act as resonant scattering centers for the phonons. 
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